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Population transfer between low lying Rydberg states independent of the initial state is realized 
using a train of half-cycle pulses with pulse durations much less than the classical orbit period. 
We demonstrate experimentally the transfer of population from initial states around n — 50 with 
10% of the population de-excited down to n < 40 as well as up to the continuum. This is the first 
demonstration of a state- independent de-excitation technique applicable to the currently produced 
state distribution of antihydrogen. The measured population transfer matches well to a model of 
the process for 1-D atoms. 



I. INTRODUCTION 

The control of atoms in Rydberg states is impor- 
tant for many different applications including the study 
of ant i- hydrogen, the study of wavepackets in Rydberg 
atoms, and the use of excitation blockades for quantum 
computing p], d, Precise tests of CPT using antihy- 
drogen requires trapped atoms in the ground state, but 
currently atoms are produced in a broad range of excited 
states. An effective means of de-excitation is thus re- 
quired to ensure atoms reach the ground state and are 
trapped before annihilating on the trap walls. Stimu- 
lated de-excitation using narrow-band lasers is clearly in- 
feasible due to the many thousands of states populated 
during three-body recombination. Hence an alternative 
state-independent technique is required. Efficient de- 
excitation of anti-hydrogen atoms compatible with the 
techniques used currently to produce anti-hydrogen re- 
quires a scheme that meets several requirements. First 
it must be capable of simultaneously de-exciting a range 
of states whose radius varies from 0.5 jam to less than 
100 nm; this range is equivalent to states in no magnetic 
field with principal quantum numbers, n, between n < 35 
and n = 90 [4]. Second the current techniques to produce 
antihydrogen atoms result in a quasi-static production 
rate over at least several seconds; a de-excitation scheme 
must thus not require a specific timing sequence. Lastly, 
atoms with thermal velocities comparable to the 1.2K 
production environment reach the walls of the trap en- 
closure in under 250 /is which sets a limit for the required 
de-excitation rate. 

Several groups have proposed schemes using trains 
of ultra-short unipolar electromagnetic pulses (half-cycle 
pulses or HCP) to accomplish this de-excitation. When 
the length of the pulse is much shorter than the orbit pe- 
riod of the electron, T or ut oc n 3 , these pulses effectively 
kick the electron transferring momentum regardless of 
the initial state [5]. Depending on the phase of the kick 
relative to the electron's orbit this momentum transfer 
can either drive the atomic state up towards the contin- 



* speck@rowland.harvard.edu 



uum or down towards the ground state. Hu and Collins 
solved the applicable time-dependent Schrodinger equa- 
tion which predicts that a fixed repetition rate train of 
pulses should result in an oscillating population distri- 
bution in number of half-cycle pulses with excursions 
from the initial state of over An = ±10 [6). Kopyciuk 
and Parzynski noted that, in a model system based on 
1-D atoms, by changing the repetition rate of the half- 
cycle pulses during an experimental sequence atoms can 
be driven either up or down in principal quantum by 
decreasing or increasing the repetition rate respectively. 
Their model however does not predict the oscillatory phe- 
nomenon seen by Hu and Collins 0,0. 

Experimentally, half-cycle pulse trains have been 
demonstrated to coherently transfer Rydberg atoms be- 
tween n-states of 350 and 700 [9]. As the orbit period 
for these atoms is on the order of tens of nanoseconds, 
electrical impulse generators are used to generate pulses 
with temporal widths less than the orbit period. However 
this experimental technique cannot be extended to atoms 
with n < 100 due to the much shorter pulse lengths re- 
quired when the orbit periods are less than 100 ps. Single 
pulses to atoms in states comparable to those produced 
for antihydrogen study have been shown to cause popu- 
lation redistribution to nearby n and / states but ioniza- 
tion occurs at high pulse amplitudes instead of inducing 
larger changes in quantum state [TO, EU- These systems 
utilize amplified femtosecond lasers with repetition rates 
limited to less than around 10 kHz to produce half-cycle 
pulses; they thus cannot be extended to produce pulse 
trains to drive large changes within the lifetime of a Ry- 
dberg atom. Utilizing adiabatic rapid passage through 
a series of microwave resonances has also been demon- 
strated as a technique to drive population transfers but 
requires specific knowledge of the dressed atomic states 
which is difficult in strong inhomogeneous magnetic fields 
such as those used to trap antihydrogen [12l 1 131] . 

In this paper we report on the first experimental 
demonstration of a technique that meets all of the re- 
quirements for efficient de-excitation. An 82 MHz rep- 
etition rate photoconductive switch produces a train of 
primarily unipolar pulses with pulse widths on the order 
of 1 ps and amplitudes greater than 500 V/cm allow- 
ing for the manipulation of atoms with n > 20 through 
impulsive kicks from these pulses. By applying thou- 
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sands of pulses within the radiative lifetime of an atom, 
T ~ 100 jis for n ~ 40, over 10% of the initial population 
in n = 52 can be driven to n = 40. As radiative lifetimes 
scale strongly with n, r ra d oc n 5 / logn, this reduction in n 
results in a factor of 4 reduction in lifetime [14]. The pop- 
ulation transfer follows the predictions of the 1-D model 
system used in the theoretical work but instead of a co- 
herent transfer from each pulse the population decoheres 
after each pulse. 

II. APPARATUS 

Fig. [TJa) shows the primary apparatus used in this ex- 
periment. Approximately 3 x 10 8 85 Rb atoms are trapped 
in a 4 mm diameter cloud in the center of a dispenser 
loaded magneto-optical trap. The combination of the 
primary cooling laser at 780 nm and a tunable 480 nm 
optical parametric oscillator producing 2 mJ, 10 ns long 
pulses with linewidth of 4 cm -1 excite 1 in 10 4 atoms 
in two steps from the 5S ground state to Rydberg nS 
or nD states. Field plates of 80% transparent stainless 
steel mesh separated by 2.2 cm can be biased produc- 
ing an electric field that field ionizes the Rydberg atoms. 
An electron multiplier detects these ionized electrons. By 
applying a linear ramped field in time the entire field ion- 
ization distribution can be mapped out in a single shot. 
This field typically ramps from to 700 V/cm in 10 /as. 

Half cycle pulses are generated within the vacuum 
system using a photoconductive switch fabricated on 
a semi-insulating GaAs wafer with micro-fabrication 
techniques [HI . SI-GaAs is a slow photoconductive switch 
with a carrier lifetime of more than 100 ps, a character- 
istic which is necessary for producing half-cycle pulses 
(HCP) as we will discuss later. The semiconductor wafers 
used here are cut along the (100) direction and have a 
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FIG. 1. Schematic of the apparatus. Atoms trapped in the 
center of the magneto-optic trap are excited to a Rydberg 
state, interact with a train of half-cycle pulses generated by 
the emitter antenna on the left side, and then are detected by 
field ionization from the field plates surrounding the central 
trapping region. The emitter and photoconductive detector 
geometries are shown in (b) and (c) respectively. 
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FIG. 2. Half-cycle pulse waveform for a biasing electric field 
of 2 kV/cm. The dashed line shows the measured current 
from the photoconductive detector while the solid line is the 
derivative of the current which is proportional to the radiated 
electric field. 

resistivity of 10 7 Q • m. Aluminum strip lines with a gap 
of 700 jam for biasing the switch are directly deposited 
on GaAs using electron-beam evaporation to a thickness 
of between 0.5 jam and 0.8 jam. Fig. DJb) shows the ba- 
sic geometry. A similar wafer is used as a detector but 
instead of a stripline antenna a dipole antenna with a 
gap of 6 jam is fabricated on the detector as shown in 
Fig. DJc). A femtosecond Tiisapphire oscillator with a 
maximum pulse width of 100 fs illuminates the emitter 
and detectors. This mode-locked laser operates around 
a center wavelength of 810 nm with a repetition rate of 
82 MHz. An average laser power as high as 450 mW fo- 
cused to a spot size of 50 /am was used to pump the 
emitter. 

To compensate for the strongly diverging nature of the 
emitted HCP radiation, a combination of a 10 mm di- 
ameter hemispherical silicon lens located on the back 
of the emitter wafer along with a pair of plano-convex 
lenses are used to refocus the pulses at the center of the 
trapped atom cloud. The two lenses have focal lengths of 
35 mm and then 50 mm and are respectively composed 
of TPX (polymethylpentene) and high-density polyethy- 
lene which are both transparent in the terahertz region. 
They are placed such that all frequencies are refocused at 
the center of the trapped atoms. An identical sequence 
of lenses directs and refocuses the outgoing pulses on 
the photoconductive detector located outside the vacuum 
system. A high-resistivity Si window allows the HCP 
pulses to exit the UHV system. 

Fig. [2] shows typical half-cycle pulses with a temporal 
FWHM of approximately 1 ps and a peak field of over 
500 V/cm produced by our setup. The fast excitation of 
carriers in the GaAs wafer due to the pump pulse cou- 
pled with the strong bias electric field results in a sudden 
current surge. Since in the dipole approximation the far 
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field radiation is proportional to the time derivative of 
the current, this current surge produces a strong positive 
pulse while the long carrier lifetime causes a slow decay 
in the current leading to a long negative tail. After prop- 
agation in free space, the strong divergence of the low 
frequency components of the pulse as well as the Gouy 
phase induced by focusing converts this pulse into the 
symmetric pulse shown in Fig. [2l 

The pulses are detected using the inverse technique in 
which the half-cycle pulse electric field is used to induce a 
current within the detector. The measured current den- 
sity is then given by j = e\i f^°EHCp(t)n(t)dt where 
e is the charge of the carriers, \i is the carrier mobility, 
and n(t) is the carrier density. For a semiconductor like 
SI-GaAs with a long carrier lifetime, n(t) can be approx- 
imated by a step function at time r, corresponding to 
the arrival of the probe pulse, so this equation reduces to 
j(r) oc EHCp(t — T)dt or equivalently the electric field 
of the half-cycle pulse is proportional to the derivative of 
the measured current. Using the simplifying assumption 
that the carrier mobility at terahertz frequencies is the 
same as that at DC frequencies and accounting for both 
the exponential decay in carrier density caused by the 
approximately 300 ps carrier lifetime of SI-GaAs and the 
reflection losses at the interface between the detector and 
free space, the current induced by a known DC biasing 
field is used to calibrate the detector. This calibration 
has been checked by comparing this detected signal with 
an electro-optic detector with a known absolute calibra- 
tion. 

Each experimental run is synchronized relative to a 
laser pulse from the OPO which provides the excitation 
to a given Rydberg state. 3 ms before the pulse the trap- 
ping quadrupole magnetic field is turned off allowing all 
resulting eddy currents to dissipate prior to the excitation 
and ensuring minimal interference from the quadrupole 
field. 5 fis prior to the pulse the photoconductive switch's 
bias field is switched on allowing it to stabilize prior to 
any Rydberg atom excitation. This eliminates any effects 
from lower amplitude pulses being emitted as the field 
is being raised to the full value. Approximately 30,000 
atoms are excited to a Rydberg state at time t = 0. Af- 
ter a variable time between and 150 /is, the half-cycle 
pulse biasing field is turned off with a fall time of 30 
ns. The final state distribution is measured by apply- 
ing a 10 /is long ionizing field ramp at a fixed time of 
150 /is after the excitation pulse. By fixing the ioniza- 
tion field ramp in time relative to the excitation pulse, 
we ensure that no effects due to radiative decay or transi- 
tions induced by blackbody radiation are convolved into 
the measured state distribution because the only change 
between experimental runs is the length of time the half- 
cycle pulses are applied. Fig. [3] (e) and (f) show the 
measured normalized distribution for a fixed number of 
half-cycle pulses as the delay between the Rydberg exci- 
tation laser pulse and the ionization ramp is varied. At 
long delay times there is a small decrease in the number 
of atoms ionized at high fields as expected due to these 
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FIG. 3. Left plots: Normalized fraction of the atomic popula- 
tion that ionizes at a specific electric field as a function of the 
number of applied half-cycle pulses, (a) shows experimental 
data and (c) shows the calculated values for the 1-D model 
with a momentum transferred per pulse of q — 7.9 x 10 -5 
in atomic units. Right plots: Atomic state distribution based 
on conversion from measured ionizing field using the expected 
fields for red-shifted Stark states, (b) shows the experimen- 
tal data while (d) is calculated using the 1-D model. Bottom 
plots: Evolution of the state distribution in time after 6500 
(e) or 2400 (f) HCP's are applied. The evolution is measured 
by delaying the ionizing ramp relative to the initial Rydberg 
excitation pulse. The dashed lines shows the ionization ramp 
delay of 150 /isused for all other plots. 

states having reduced lifetimes but at a delay of 150 /is, 
used for all of our other data, there is a negligible effect 
on the overall state distribution. 



III. RESULTS AND DISCUSSION 

Fig. [3ja) shows the normalized state distribution mea- 
sured as a function of the number of half-cycle pulses. 
Initially all atoms ionize near 80 V/cm as set by the ex- 
citation wavelength and resulting state. As impulsive 
kicks are applied to this initial distribution the popula- 
tion spreads out in a diffusive manner. During this pro- 
cess some atoms are lost due to field ionization by the 
small static lV/cm clearing field or from radiative decay 
to states that are not ionized by the maximum field of 
700 V/cm applied during the detection ramp. To map 
between ionizing electric field and the approximate prin- 
cipal quantum number distribution shown in Fig. [3](c) we 
utilize the electric field at with the electron is no longer 
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FIG. 4. (a) Field ionization trace showing peaks from popula- 
tions in many different principal quantum numbers. Dashed 
lines are the calculated ionizing field for red Stark states, (b) 
Ionization peak locations versus principal quantum number. 
Note the close agreement to the red Stark state calculated 
value of F = (l/9)/n 4 . The dashed and dotted lines shows 
the expected ionization fields for nS and nD states, respec- 
tively, following adiabatic pathways towards ionization, (c) 
Similar fit assuming adiabatic ionization of nD states. 



classically bound. In atomic units this field is F = 
Here A = 1/16 for adiabatic ionization or A = 1/9 for 
red-shifted Stark states that diabatically ionize which we 
employ for data analysis from now on[16]. Justification 
for the choice of red-shifted Stark states is shown in Fig. 
IU Part (a) shows a field ionization trace after 1600 half- 
cycle pulses have been applied to the population. More 
than 10 peaks corresponding to discrete atomic states are 
clearly visible. Part (b) plots these peak positions versus 
the best fit principal quantum number and corresponds 
closely to a distribution from red-shifted Stark states. 
The inset shows a comparison to the expected distribu- 
tion if we assume nD states have been produced which 
does not agree as well. The predominant excitation of 
red-shifted Stark states is due to the small electric field 
present during excitation[17j. Due to the large quantum 
defects of I < 2 states in Rb, these optically accessible 
states are displaced at zero field from the Stark manifold 
of states with I > 2. While the first Stark manifold that 
a nS state encounters is the n — 3 manifold with which 
there is only a small coupling due to the large difference 



in n, a nD state encounters the downhill, red-shifted state 
of the — 1 Stark manifold where there is a much larger 
coupling. At this point much of the D character is spread 
to these red-shifted states and it becomes possible to op- 
tically excite them when the electric field is high enough 
to have reached the first crossing, F ~ l/3n 5 . 

Assuming the length of the half-cycle pulse is short rel- 
ative to the electron orbit time, the transition amplitude 
from initial state, z, to final state, /, due to an interaction 
with a HCP can be reduced to the inelastic form factor 
of the atom, Tf^ = (f\e~ zq ' r \i) where q is the momentum 
transferred to the atom. With an initial density matrix, 
p(0), given by a statistical mixture of the measured initial 
state distribution, the state distribution resulting after j 
pulses is given by 



P(j) = (e- iHt <-T) j p(0) (Th 



iHt r 
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-l/2n 2 



where t r is the period of the pulse train and H 
is the unperturbed Hamiltonian of the system. 

The applicable form factor has been analytically cal- 
culated for hydrogenic Stark states with parabolic quan- 
tum numbers rii, 77,2, m, and n = n\ + ri2 + |m| + 1 
where the atomic quantization axis and the HCP electric 
field polarization are parallel[18]. In this configuration 
the only applicable selection rule is Am = 0. When the 
quantization axis and the HCP electric field polarization 
are perpendicular, there is currently no simple analyt- 
ical formula. In the case of weak momentum transfer, 
Tf t i ~ Si j — iq x (f\x\i) and the applicable matrix ele- 
ments are known [19j. In contrast to the parallel case, 
the selection rule for the magnetic quantum number is 
now Am = ±1. As red-shifted Stark states are predomi- 
nantly one-dimensional atoms with the expectation value 
of the axial electron position much greater than the ra- 
dial position, the interaction can be further simplified to 
that of 1-D atoms interacting with an impulsive kick of 
given momentum transfer, q. The transition amplitude 
for a single impulsive kick in the 1-D approximation is 
given by [20]: 
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where n and n 1 are the initial and final principal quantum 
numbers, F(a, 6, c, y) is the hypergeometric function, A = 
n _1 + (n / ) _1 — iq, and z = —An'n [(n — n') 2 + (qn'n) 2 ] 1 . 

Fig. [5] compares the probability, |T/^| 2 , of a tran- 
sition for these three different cases from a red-shifted 
Stark state where n\ = and |m| < 5 to a final state 
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FIG. 5. Transition probability from n — 62 to n — 62 + An 
for a momentum transfer of q = 7.9 x 10 -5 . The probabil- 
ities from the 1-D model (Eq. [2} and from 3-D red-shifted 
Stark states with \m\ < 5 kicked by a pulse with electric field 
polarization along the quantization axis are almost the same. 
The probabilities from the same Stark states kicked instead 
by a pulse with electric field polarization perpendicular to the 
quantization axis are more than an order of magnitude lower. 



with a given change in principal quantum number. Two 
points from this figure are important. First the over- 
all shape of the curve is similar for all three cases but 
the probability of a transition in the perpendicular kick 
case is approximately an order of magnitude lower. Since 
AE = (2po • q + q 2 )/2 where po is the initial momen- 
tum and with perpendicular polarization po • q = 0, one 
would expect the reduction in energy transfer and thus 
transition amplitude when kicked in this configuration. 
Second, in the two 3-D cases the transition amplitude is 
clustered around similar states. Hence the likelihood of 
driving a transition away from red-shifted Stark states is 
many orders of magnitude lower. Given these results, the 
1-D model is expected to apply quite closely apart from 
the correction in overall amplitude necessary for the per- 
pendicular kick case as atoms that start in a primarily 
1-D state will stay in similar states. 

Unless a fast decoherence term is added, this model 
shows interference effects and emphasizes eigenstates 
which have stationary phase at the frequency of the pulse 
train. This is not seen in the actual experimental data 
suggesting that decoherence occurs on a time scale simi- 
lar to the time between half-cycle pulses. This is a shorter 
time scale than dephasing due to blackbody transitions 
or spontaneous emission. Also no dependence on density 
is observed precluding a collisional process being respon- 
sible for the decoherence. The most probable reason for 
the fast decoherence is the fact that the polarization of 
the half-cycle pulse is nearly perpendicular to the quan- 
tization axis set by the static external electric field. Each 
kick must change m by ±1 allowing for many more states 
to be populated and small changes in the time between 
HCP's will result in dephasing between these states. 

Fig. [3](d) shows the results of applying the above 1-D 
model to the initial state distribution shown in Fig. [3](b) 
measured when no half-cycle pulses have been applied. 
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FIG. 6. Electric field for which the given percentage of the 
atomic ensemble has ionized as a function of the number of 
half- cycle pulses applied. The initial state of the atoms varies 
from (a) n — 62, (b) n — 52, to (c) n — 46. Solid lines 
show the expected value from the 1-D model with momentum 
transfer per pulse of q = 7.9 x 10 -5 in atomic units. 



The predicted ionizing field distribution shown in Fig. 
(He) is calculated by assigning an electric field distribu- 
tion to each principal quantum number centered at the 
expected field for red-shifted Stark states with a width 
proportional to the distance between adjacent state field 
ionization center values. 

To compare the model and experimental data in more 
detail, Fig. [6] shows the electric fields at which 10%, 
50%, and 90% of the atoms have ionized for a series 
of different initial states. The agreement between the 
experimental data and the one free parameter model is 
quite good especially for atoms in higher n states. The 
lack of agreement for lower n states, or equivalent ly pop- 
ulations that ionize at higher electric fields, in Fig. [6] 
is caused by a series of factors. First at lower princi- 
pal quantum numbers the 1-D approximation, even for a 
red-shifted Stark state, begins to break down. Addition- 
ally half-cycle pulses applied to atoms in these states are 
more likely to drive atoms to states for which the 1-D 
approximation no longer applies. In addition the width 
in electric field over which a state ionizes is proportional 
to the field at which it ionizes. Thus for lower n states 
the width is much greater resulting in a lower signal at 
a given field and increasing the sensitivity to noise and 
the exact details of how the conversion between state and 
ionizing field distributions is performed. 

The maximum momentum transferred per pulse uti- 
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FIG. 7. Best fit for the impulse in atomic units transferred by 
a single half-cycle pulse to the 1-D model described in the text, 
(a) the initial state of the atoms was held fixed at n — 62 while 
the amplitude of the half-cycle pulse was scanned. The solid 
line shows the measured amplitude of the half-cycle pulse, (b) 
the half cycle amplitude was held constant at a biasing field 
of 4.3 kV/cm while the initial atomic state was varied. 



lized in the model was calculated by a best fit to the 
data for the field at which 10% of the total population 
was ionized including several initial state distributions 
as shown in Fig. [T^b). Averaging over these points, 
we find a calculated momentum transfer per pulse of 
q = 7.9 x 10 -5 . From Eq. [2] and this calculated mo- 
mentum transfer, there is less than a 1% probability of 
an atom changing principal quantum number by ±1 from 
the action of a single pulse for states near n ~ 50 and a 
negligible probability to change by higher amounts. Con- 
sequently this model behaves as a diffusive process and 
the spreading of population density is expected. From 
the measured electric field of the half-cycle pulse we ex- 
pect a maximum momentum transferred by a single pulse 
of 5 x 10 -4 in atomic units. This is a factor of 5 times 
more than the maximum calculated from the 1-D model 
fit as shown in Fig. [7(b). However, correcting the calcu- 
lated value for the reduced transition amplitudes when 
the HCP kicks the atom perpendicular to its quantization 
axis, we find <2~3xl0 -4 which is in close agreement. 



IV. CONCLUSION 

In summary, a train of half-cycle pulses produces large 
diffusive population transfers independent of the initial 
n-state of the Rydberg atoms and operates regardless of 
the time structure creating Rydberg atoms. These trains 
are thus ideal for the control and efficient de-excitation of 
Rydberg antihydrogen atoms as they meet all of the con- 
ditions required for compatibility with antihydrogen pro- 
duction discussed previously. The ability to accurately 
model the population transfer demonstrated in this text 
utilizing analytical calculations of the inelastic form fac- 
tor will allow for the efficient development of enhanced 
protocols utilizing, for example, chirped pulse trains to 
both enhance the efficiency of the transfer as well as to 
drive atoms toward a desired final state. These protocols 
will depend on coherent population transfers and thus 
the rapid dephasing observed in our experiments must 
be reduced prior to their implementation. 

Regardless of further developments, the population 
transfers demonstrated in this paper will greatly increase 
the probability of trapping antihydrogen atoms. The 
10% increase in /-mixed population with n < 40 has 
a radiative lifetime a factor of 4 shorter than that of 
n = 52; furthermore, during the 150 /is that half-cycle 
pulses were applied, an additional 2% of the population 
reached the ground state. For antihydrogen applications, 
this reduction from over a 4 ms to a 1 ms lifetime should 
increase the population that reaches the ground state 
before annihilating on a trap wall by over an order of 
magnitude based on the current ATRAP experimental 
configuration [2lj . 
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